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Abstract 

The cumulants and factorial moments of photon distribution for 
squeezed and correlated light are calculated in terms of Chebyshev, 
Legendre and Laguerre polynomials. The phenomenon of strong os- 
cillations of the ratio of the cumulant to factorial moment is found. 
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1 INTRODUCTION 

The coherent light [1],[2],[3] has the photon distribution function described by the standard 
Poisson distribution. The nonclassical states of light, for example squeezed states [4], [5], [6], 
Schrodinger cat states [7], [8], [9], correlated states [10] have the photon distribution functions 
which differ essentially from poissonian ones demonstrating either superpoissonian or sub- 
poissonian behaviour which distinguishes the nonclassical types of light from the coherent 
light considered by definition as classical one. 

One of the important differences of photon distribution function of the nonclassical types 
of light from the poissonian distribution of the coherent light is the possible strong oscillations 
of the photon distribution function. For the squeezed light these oscillations have been found 
in [11], [12]. For the correlated light the existence of such oscillations has been demonstrated 
in [13]. For even and odd coherent states (Schrodinger cat states) the fast oscillations of the 
photon distribution function are connected with the absence of the states with even numbers 
of photons in odd cat states and the states with odd numbers of photons in even cat states 
[7], [14], [15]. Such properties of the photon distribution oscillations for the nonclassical types 
of light are preserved also for the multimode electromagnetic radiation. It was shown for 
two-mode squeezed light in [16], [17], [18] and for even and odd cat states in [19], [20]. 

On the other hand, in high energy physics another characteristics of the multiparticle 
distribution function has been widely used, namely, factorial moments and cumulants (for 
the review see [21]). It was shown that the cumulants and the functions of these quantities 
are very sensitive to the details of the particle distribution demonstrating the oscillating 
behaviour [22], [23], [24]. The properties of the cumulants of the integer rank for the known 
distributions as poissonian and binomial ones are well known and for these distributions the 
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oscillatory behaviour is absent [24]. 

The behaviour of the cumulants (in particular, their oscillations) for the squeezed states, 
correlated states or other nonclassical states has not been studied till now in all details. At 
the same time, for the most general photon distribution function corresponding to squeezed 
and correlated one-mode light at finite temperature the explicit expression in terms of Her- 
mite polynomials of two variables has been obtained in [25] and it was widely used in 
[26], [27], [28], [29]. 

The aim of this work is to obtain the explicit expressions for the cumulants and factorial 
moments of the photon distribution function for the squeezed and correlated light at finite 
temperature. We demonstrate that the cumulants possess the strongly oscillating behaviour 
in the region of slight squeezing where the photon distribution function itself has no oscilla- 
tions. And vice versa in the region of large squeezing, where the photon distribution function 
strongly oscillates, the cumulants behave smoothly. Thus the behaviour of the cumulants may 
provide very sensitive method of detecting very small squeezing and correlation phenomena 
due to presence of strong oscillations. 

2 PHOTON DISTRIBUTION FUNCTION 
AND ITS MOMENTS 

Since in what follows it will be necessary to use some characteristics of one-mode squeezed 
light, we will briefly review the main results obtained (see, for example [26]). 

Let us consider the most general mixed squeezed state of one-mode light described by the 
Wigner function W(p, x) of the generic Gaussian form with five real parameters, (x), (p), 
&xx, Vpp, c P x (first two parameters are means of position and momentum, others are matrix 
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elements of the dispersion matrix for the position and momentum), 



W(p,x) = tT^exp {-(2d)- 1 [a xx (p - (p)) 2 
+ a pp (x - (x)) 2 - 2a px (p - (p))(x - (x)j\ } , 



where 



7 2 



is determinant of the dispersion matrix. 

For the photon distribution function the following formula was obtained [25], [26] 
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is the probability to have no photons, 



= O m + (7; 
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is the trace of the dispersion matrix, H^ n - Hermite polynomials of two variables. 
Elements of the symmetric matrix 



R = 



Rn R12 
R\2 R22 



determining the Hermite polynomials are given by the formulas: 



and two arguments of the polynomials are defined by the equation: 

2 [(T - l)z* + (<7pp - a xx + 2ia px )z] 

y-i — v — — . 
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The complex parameter z is given by the relation: 

z = 2- l ' 2 ((x)+i(p)) . 
The generating function for the photon distribution function was also obtained in 
G(u) = P 
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The photon distribution function is related to as follows: 



n\ du n 



u=0 



The above-mentioned normalized cumulants and factorial moments are defined as 

1 cMnG(w) 
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respectively, and related by the following recursion relation: 
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where 
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are the binomial coefficients. Some formulas for the cumulants of photon distribution function 
can be found in [30], [31]. The factorial moments have been widely used in particle physics [21] 
to analyze intermittency properties [32] of fluctuations. However, they are less instructive 
than the cumulants or their functions. It was shown in [22], [24] that the ratio of cumulant 
to factorial moments, i.e. the function 



is a very sensitive measure of tiny details of the multiplicity distribution. In particular, it 
can be used to distinguish between different distributions which otherwise look quite similar. 



It was already mentioned in the Introduction that the photon distribution function exhibits 
an oscillatory behaviour if we deal with highly squeezed states ( T = a pp +a xx 1 ) for large 
values of the parameter z (2.3). A question arises: is it possible to obtain a similar "abnormal" 
behaviour of other characteristics of the photon distribution, namely, cumulants, factorial 
moments and the function H q defined in (2.8)? If yes, then in what region of parameters of 
the function (2.1) such anomalies can take place? The present section is dedicated to the 
solution of this problem. 

The direct differentiation of the function \nG(u) at u — 1 yields: 



H q = KJF, 



(2.8) 



3 ANALYSIS OF K q , F q AND H q 
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where [26] 
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High oscillation are obtained at strong squeezing (large T) or at large values of \z\ 2 . 
It is known [27] that 

T > 1 , d > \ . (3.2) 
Using (3.2) the following inequalities can be easily obtained: 

Ai > 1 and 

A 2 < or A 2 > 1 . (3.3) 

Expression in the square brackets in (3.1) consists of two terms: 

1 f 1 , CiAi \ 



(Ax - l) 9 V 2 1 - A i 



(3.4) 
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The first term is of constant sign. The second one is oscillating in the case A 2 < 0. With aim 
to obtain the oscillations of the whole function K q we will treat only this case: 

A 2 < . 

Then 

1 (-!)« 



(A 2 -l)9 (1 + |A 2 |)« 

However, then it follows that 

|A 2 | > Ai > 1 , 



and the alternating term diminishes faster than the constant sign term. The terms 1 / 1 A2 — 1 j 
and l/(Ai — 1) are most close to one another if 

d — - (the pure state), 

that is used in the following. 

First of all we consider the simplest case when the value of z as given by (2.3) equals to 
zero. Then 

where 

a 7 1 T T-l 
p = d + - — — - , a 
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T q (a) and P q (a) are the Chebyshev polynomials of the first kind and Legendre polynomials, 
respectively. Let us note that the arguments of polynomials are purely imaginary but the 
whole expressions for moments are real, surely. 
For H q we obtain the expression: 



Taking d — 1/4 we have: 



H q = ■ (3-7) 
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In this case the curve H q has step-like shape at (T — 1) — > 0; steps become smoothed as 
T grows (fig.l). We should note that direct limit T — > 1 shows the discontinuous character 



of the function H q (T) at T = I. The point is that at T = 1 we are dealing with the usual 
Poisson distribution (let us remind that we treat a case d = 1/4, \z\ = 0), where H q = 5 q ±, i.e. 
Hi — 1, H q — at q ^ 1, which differs from the behaviour of if 9 at (T — 1) = 10~ 5 , depicted 
in the fig 1. The particular values of the second rank moments are very high what reveals 
extremely wide distribution (so wide distributions are unknown in particle multiproduction, 
for example). Namely, one can easily show that F 2 > 3, K 2 > 2, H 2 > 2/3. 

Consider now the case \z\ ^ 0. Since the photon distribution function is invariant with 
respect to rotation in a phase space, without loss of generality we can consider a xx = a pp 
(vpx 7^ - correlated state). By appropriate choice of the phase of (2.3) ((x) = — (p)) we 
cancel the linearly growing term q£i\±/ (1 — Ai) in (3.4). Moreover, the analogous linear term 
^2^2/ — A2) in (3.5) becomes maximal at fixed \z\. Thus we have left only two variable 
parameters T and \z\, and formula (3.1) has the following final form: 
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where 
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In the case of large T (highly squeezed state) we can obtain the finite number of oscillations 
of K q taking large value of \z\. However, the average number of photons in corresponding 
states is large, and the amplitude of the oscillations decreases exponentially due to the factor 
l/(n) q . Remind that in this very case the strong oscillations of the photon distribution 
function can be observed. 

Now let us consider the case of the slightly squeezed state, y — (T — 1) <^ 1, when photon 



distribution function does not oscillate. Impose also an additional condition 



|2 



\ z \ 

that makes possible to obtain approximate formulas for the functions K q , F q and H q . For 
K q we have the following approximate expression: 

K q = g!(-l)«-y-« . (3.9) 

Then recursion relation (2.7) yields: 

F q = q\(-iy 1 - q L q \ 1 ) , (3.10) 

where L~ l (x) are generalized Laguerre polynomials. For H q with g C7 we have: 

H q = K q /F q = -jp^fi « (-i) q+1 qh 1 - q « 1 • (3.H) 

(If 7 ^> q, the term with the highest power of 7 dominates over the rest of the sum in L" 1 ^), 
and F q — > 1 as for Poisson distribution). The exact shape of the function H q is shown in 
fig. 2. The distribution function P n does not oscillate (fig.2 a ). 

However, the most abrupt oscillations of the functions K q and H q have been obtained when 
(T — 1) <C 1, but condition 7 ^> 1 is not valid. The corresponding curves are shown in the 
figs. 3, 3 a . Note that the photon distribution function is smooth again being approximately 
equal to zero at q ^ 1. 

The most regular oscillating patterns of K q and H q are seen at (T — 1) ~ 0.1, \z\ ~ 1 
(figs.4, 4 a ). 

The alternating sign cumulants are typical also for the fixed multiplicity distribution, i.e., 
for P n = 5 nno (no = const) [24] . Let us note that there exist smooth multiplicity distributions 
which give rise to cumulants oscillating with larger period (see [22]). 



Finally we consider the opposite case when the photon distribution function P n exhibits 
strong oscillations while K q and H q behave smoothly. Such a behaviour is typical at T ~ 100, 
\z\ ~ 1 when K q exponentially grows while H q monotonically decreases with q (fig. 5). 

4 CONCLUSION 

We have shown that cumulants of the photon distribution function and their ratio to the 
corresponding factorial moments exhibit oscillating behaviour in the case of slightly squeezed 
states: (T — 1) = (a pp + a qq ) — 1 1. We have considered also the pure state with d — 1/4. 
Oscillations of the photon distribution function are absent in that case. 

We are yet unable to establish direct correspondence between the oscillations of cumulants 
and the behaviour of the photon distribution function. Somehow it should depend on the 
range of parameters considered. For example, it is known [26] that parameter d = cr pp a xx —a px 
characterizes the temperature of the system and the temperature equals to zero at d — 1/4. 
Is it possible to observe oscillations of cumulants and of the function H q at other parameters 
of the function (2.1), for example, at d > 1/4? This problem has not been studied in all 
details yet. However, evidently, the oscillations must be smoothed as d grows since it is 
analogous to the behaviour of the function P n , whose oscillations disappear with increase of 
temperature. The extcusion of the above procedure to non-integer ranks is straightforward. 
We hope to consider these problems in subsequent publications. 
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FIGURE CAPTIONS 

Figure 1: The behaviour of the function H q defined in (2.8) at d — 1/4, \z\ = 0; parameter 
T is varied: (1) (T — 1) = KT 5 , (2) T = 1.1, (3) T = 1.2. 

Figure 2: The behaviour of the function H q at d = 1/4, T = 1.01; the curves (1) and (2) 
correspond to the values \z\ 2 = 1.01, 0.8, respectively. 

Figure 2 a : The photon distribution function at d — 1/4, T — 1.01; the curves in order of 
the lowering maxima correspond respectively to the values \z\ 2 = 1.01, 0.8. 

Figure 3: The cumulants of the photon distribution function at d — 1/4, (T — 1) = 10~ 5 , 
\z\ 2 = 0.01. 

Figure 3 a : The behaviour of the function H q at d = 1/4, (T — 1) = 10~ 5 ; parameter \z\ is 
varied: (1) \z\ 2 = 0.01, (2) \z\ 2 = 0.005, (3) \z\ 2 = 0.01. 

Figure 4: The cumulants of the photon distribution function at d — 1/4, T — 1.1, \z\ 2 = 
1.1. 

Figure 4 a : The behaviour of the function H q at d = 1/4, T = 1.1; parameter \z\ is varied: 
(1) \z\ 2 = 2, (2) \z\ 2 = 1.1. 

Figure A b : The photon distribution function at d = 1/4, T = 1.1; the curves in order of 
the lowering maxima correspond to the values \z\ 2 = 2, 1.1, respectively. 

Figure 5: The smooth curve for the function H q and the oscillating photon distribution 
function P n at d = 1/4, T = 100, \z\ = 1. 
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